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Abstract 

We present an analysis of hadronic spectroscopy for Wilson valence quarks 
with dynamical staggered fermions at lattice coupling 6/g 2 = {3 = 5.6 at sea 
quark mass am q = 0.01 and 0.025, and of Wilson valence quarks in quenched 
approximation at (3 = 5.85 and 5.95, both on 16 3 x 32 lattices. We make 
comparisons with our previous results with dynamical staggered fermions at the 
same parameter values but on 16 4 lattices doubled in the temporal direction. 
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1 Introduction 



Calculations of hadron spectroscopy remain an important part of nonperturbative 
studies of QCD using lattice methods. (For a review of recent progress in this field, 
see Ref. |].) We have been engaged in an extended program of calculation of the 
masses and other parameters of the light hadrons in simulations which include the 
effects of two flavors of light dynamical quarks. These quarks are realized on the lattice 
as staggered fermions. We have carried out simulations with lattice valence quarks in 
both the staggered and Wilson formulations. Our reasons for performing simulations 
with Wilson valence quarks is twofold: First, we are interested in seeing if there are 
any effects of sea quarks on the spectroscopy of systems containing either realization 
of valence quark. Thus this work complements our parallel studies of spectroscopy 
with staggered valence and sea quarks, and of spectroscopy with Wilson valence and 
sea quarks. Second, we are interested in exploring the effects of sea quarks on simple 
matrix elements such as the pseudoscalar meson decay constant. Most previous work 
has been done with Wilson valence quarks in quenched approximation. We consider 
that mixing the two realizations is not inappropriate for a first round of numerical 
simulations. 

These simulations are performed on 16 3 x32 lattices at lattice coupling 6/g 2 = (3 = 
5.6 with two masses of dynamical staggered fermions, am q = 0.025 and am q = 0.01. 
These are the same parameter values as we used in our first round of simulations 
However, the first set of simulations have two known inadequacies. The first is that 
most of our runs were carried out on lattices of spatial size 12 3 . A short run on 16 3 
lattices with dynamical quark mass 0.01 showed that these lattices were too small: 
baryon masses fell by about fifteen per cent on the larger lattice compared to the 
smaller one. Thus the am q = 0.025 results from Ref. f2| are suspect and need to be 
redone. We also felt that we needed more statistics on the am q = 0.01 system. 

Second, nearly all of our running was done on lattices of size 12 4 or 16 4 ; these 
lattices were doubled in the temporal direction to 12 3 x 24 or 16 3 x 32 for spectroscopy 
studies. Doubling the lattice introduced strange structure in the propagators of some 
of the particles: the pion mass, in particular, showed strange oscillatory behavior as a 
function of position on the lattice. This behavior is almost certainly due to doubling 
the lattice and the best way to avoid this problem is to begin with a larger lattice 
in the temporal direction. 

Finally, it is an open question how much sea quarks affect the hadronic spectrum. 
In order to address this question, we have performed a set of simulations in the 
quenched approximation at lattice couplings (3 = 5.85 and 5.95, also on 16 3 x 32 
lattices. As the reader will see, our quenched results are rather similar to our results 
with dynamic fermions; apparently at the particular values of sea quark mass and 
lattice coupling where our simulations were performed, the effects of sea quarks can 
be absorbed into renormalizations of the lattice coupling and valence quark mass. 

Some of the results described here have been presented in preliminary form in 
Ref. ||. Several other papers which we are preparing for publication complement 
the results presented here: we are also preparing a paper on spectroscopy results 
with valence staggered quarks, an analysis of simple matrix elements with Wilson 
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valence quarks, a study of valence quark Coulomb gauge wave functions, and a study 
of glueballs and topology in the presence of dynamical staggered quarks. 

2 The simulations 

Our simulations were performed on the Connection Machine CM-2 located at the 
Supercomputer Computations Research Institute at Florida State University 

We carried out simulations with two flavors of dynamical staggered quarks using 
the Hybrid Molecular Dynamics algorithm || The lattice size is 16 3 x 32 sites and 
the lattice coupling (3 = 5.6. The dynamical quark mass is am q = 0.01 and 0.025. 
The total simulation length was 2000 simulation time units (with the normalization 
of Ref. J2|) at each quark mass value, after thermalization. The am q = 0.01 run 
started from an equilibrated 16 4 lattice of our previous runs on the ETA-10, that 
was doubled in the time direction and then re-equilibrated for 150 trajectories. The 
am q = 0.025 run was started from the last configuration of the smaller mass run, and 
then thermalized for 300 trajectories. We recorded lattices for the reconstruction of 
spectroscopy every 20 HMD time units, for a total of 100 lattices at each mass value. 

We computed spectroscopy with staggered sea quarks at five values of the Wilson 
quark hopping parameter: re = 0.1600, 0.1585, 0.1565, 0.1525, 0.1410, and 0.1320. 
The first three values are rather light quarks (the pseudoscalar mass in lattice units 
ranges from about 0.25 to 0.45) and the other three values correspond to heavy quarks 
(pseudoscalar masses of 0.65 to 1.5). The first three values are the ones we used in 
our first round of experiments. We computed masses of mesons with all possible 
combinations of quark and antiquark mass; we only computed the masses of baryons 
made of degenerate mass quarks. 

The quenched simulations were performed at lattice couplings of (3 = 5.85 and 
5.95, also on 16 3 x 32 lattices. They used the Hybrid Monte Carlo algorithm^ These 
simulations had a total length of 3200 time units at (3 = 5.85 and 3800 time units at 
f3 = 5.95. We recorded lattices for Wilson valence spectroscopy every 40 time units, 
for a sample of about 90 lattices at each coupling. Our quenched spectroscopy was 
done at hopping parameter values designed to reproduce as accurately as possible our 
earlier (3 = 5.6 running: we used re = 0.1585 and 0.1600 at (3 — 5.85 and re = 0.1554 
and 0.1567 at (3 = 5.95. We only computed spectroscopy for hadrons with degenerate 
quarks. 

For the spectroscopy we used periodic boundary conditions in all four directions 
of the lattice. Our lattices are long enough in the time direction and our interpolating 
fields are good enough that we are always able to extract an asymptotic mass. The 
use of open boundary conditions introduces edge effects which are hard to quantify 
and we have chosen to avoid them in the current round of simulations. 

We calculated hadron propagators in the following way: We fix gauge in each con- 
figuration in the data set to lattice Coulomb gauge using an overrelaxation algorithm 
and use sources for the quarks which spread out in space uniformly over the simu- 
lation volume and restricted to a single time slice ( "wall" sources 0). This source 
is nonzero only on sites which form one checkerboard of the lattice (the sum of x 
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plus y plus z coordinates is an even number). Our inversion technique is conjugate 
gradient with preconditioning via ILU decomposition by checkerboards]^]. We used 
a fast matrix inverter written in CMIS (Connection Machine Instruction Set) ||10|| ■ 

We use both a spread out sink as well as a "pointlike" sink where all the quarks 
lines end on the same site. The "wall" sink is identical to the source, but sums over 
all sites. We label hadron propagators with wall sources and point sinks as "WP" and 
those with a wall source and a wall sink as "WW" . We combine the quark propagators 
into hadron propagators in an entirely conventional manner. For Wilson hadrons we 
use relativistic wave functions ||1 1||. The baryon wave functions are: 



Proton: 

\P,s) = (uCj 5 d)u s 

= (u\d 2 — u 2 di + u^di — u±dz)u s 

Delta: 

|A x ,s) = (uid 2 + u 2 di + u 3 d 4 + u 4 d 3 )u s 

\A 2 ,s) = (uidz - u 2 d A + u 3 d 1 - u A d 2 )u s (1) 

We have measured meson correlation functions using spin structures ^75^ and ^7075^ 
for the pion and '073^ and ^7073'?/' f° r the rho, which we refer to as "kind = 1" and 
"kind = 2" for the pseudoscalar and vector, respectively. 

To extract masses from the hadron propagators, we must average the propagators 
over the ensemble of gauge configurations, estimate the covariance matrix and use 
a fitting routine to get an estimate of the model parameters. The lattices used for 
Wilson spectroscopy with staggered sea quarks are separated by 20 HMD time units 
and do not show any discernable time correlations with each other. The quenched 
simulations, spaced 40 Hybrid Monte Carlo time units apart, show some residual time 
correlation when we compare the error on the pion effective mass blocking various 
numbers of successive lattices together. We attempt to take these correlations into 
effect by blocking three successive lattices together before fitting the data. 

We use the full covariance matrix in fitting the propagators in order to get a mean- 



ingful estimate of the goodness of fit. Reference |12|1 discusses this fitting procedure 
in detail. 



3 Spectroscopy Results 

We determined hadron masses by fitting our data under the assumption that there 
was a single particle in each channel. This corresponds to fitting for one decaying 
exponential and its periodic partner. We calculated effective masses by fitting two 
successive distances, and also made fits to the propagators over larger distance ranges. 

In selecting the distance range to be used in the fitting, we have tried to be 
systematic. We somewhat arbitrarily choose the best fitting range as the range which 
maximizes the confidence level of the fit (to emphasize good fits) times the number 
of degrees of freedom (to emphasize fits over big distance ranges) divided by the 
statistical error on the mass (to emphasize fits with small errors). We typically 
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restrict this selection to fits beginning no more than 11 or 12 timeslices from the 
origin. 

3.1 Simulations with sea quarks 

We computed spectroscopy for the six possible values of valence quark hopping pa- 
rameter given above, with both 'WP' and 'WW correlation functions. We first show 
the global results to spectroscopy, giving the best-fit value for the mass for each value 
of hopping parameters. We display masses as a function of the average hopping pa- 
rameter + k%) for mesons and as a function of k for baryons (recall that for each 
sea quark mass we only studied baryons in which all three quarks have the same mass) 
in Figs. In all these figures masses are quoted in lattice units. We display plots of 
effective mass in Fig. § and of mass versus D min (with D max = 16) for k = .1600 data 
in Fig. Since the "kind=l" and "kind=2" operators produce essentially identical 
spectroscopy we only show "kind=l" results in these figures in an attempt to avoid 
clutter. Best fit values for each particle are shown in Tables |T]-|8|. 

For all but two cases the fitting was straightforward. However, we had two difficult 
data sets, k = .1320 spectroscopy and the am q = 0.01, k = .1600 delta. 

For all particles except those containing one or more of the heaviest quarks both 
kinds of correlation functions give consistent results, with the 'WP' correlators typ- 
ically showing the smallest uncertainty. For mesons or baryons containing a heavy 
k = .1320 quark, however, the 'WP' fits do not settle down to an asymptotic value. 
The effective mass drifts continuously with t value and fits to a range have unaccept- 
ably high chi-squared's. The WW fits are more acceptable (have a chi-squared near 
one per degree of freedom). Possibly what is happening is this: For these states the 
wall source has poor overlap on the lightest hadron in a channel, since bound states 
of heavy quarks have small spatial extent. The WP correlators do not give a varia- 
tional bound and it happens that they approach an asymptotic mass from below. The 
WW correlators approach an asymptotic value from above, but are noisier; one gets 
a statistically more acceptable fit because of larger uncertainties on the individual 
points. As an example, we show effective masses as a function of t for the ma = 0.01, 
k = 0.1320 data in Fig. 0. 

Next, the am = 0.01, k = .1600 delta mass is considerably lighter than we saw in 
our old running and nearly degenerate with the nucleon. At this light valence quark 
mass the WW operators are so noisy that they are useless, but the same behavior 
is seen in 'WP' delta operators with both spin structures of Eqn. [I]. A comparison of 
the two operators is shown in Fig. |5]d (effective masses) and |p (fits to a range d m i n 
to 16). The signal for a light delta appears to be stable for fit distances ranging from 
dmin — 4 out to 8 or 9, and then the signal deteriorates so rapidly that we cannot 
trust our fits. We do not know of a reason for this effect. A coding error would 
mix some component of the nucleon into the delta and the asymptotic mass in the 
delta channel would be the nucleon's. However, the am q = 0.025 delta is measurably 
heavier than the nucleon, which argues against a coding error. One should note that 
the WP' correlators are not variational. It is possible that we are seeing a signal 
attempting to approach an asymptotic value from below, which becomes lost in the 
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Figure 1: Best fit masses (from fits to a range) for the pseudoscalar as a function of 
the average hopping parameter. Data are labelled by type (WP or WW) (described 
in the text) by crosses (WP) and diamonds (WW). Figure (a) is for sea quark mass 
am q = 0.01, (b) for am q = 0.025. 
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Figure 2: Best fit masses (from fits to a range) for the vector meson as a function of 
the average hopping parameter. Figure (a) is for sea quark mass am q = 0.01, (b) for 
am q = 0.025. 
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Figure 3: Best fit masses (from fits to a range) for the proton as a function of hopping 
parameter. Figure (a) is for sea quark mass am q = 0.01, (b) for am q = 0.025. 



9 





0.13 0.14 0.15 0.16 




0.13 0.14 0.15 0.16 

/c 

Figure 4: Best fit masses (from fits to a range) for the delta as a function of hopping 
parameter. Figure (a) is for sea quark mass am q = 0.01, (b) for am q = 0.025. 
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Figure 5: Effective mass fits to k — 0.1600 data: (a) pion, (b) rho, (c) proton, and 
(d) delta. Data are labelled by type (WP or WW) by crosses (WP) and diamonds 
(WW). 
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Figure 6: Fits from t = D min to 16 to k — 0.1600 data: (a) pion, (b) rho, (c) proton, 
and (d) delta. Correlator types are labelled as in Fig. |l|-|5|. 
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Figure 7: Effective mass fits to k — 0.1320 data: (a) pion, (b) rho, (c) proton, and 
(d) delta. Data are labelled by type (WP or WW) by crosses (WP) and diamonds 
(WW). 
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noise before its asymptotic value is reached. 

Our data can be compared with our previous 16 4 runs: The k = .1585 and .1600 
pseudoscalar vector, and proton masses are consistent with the old numbers. The 
k = .1585 deltas are consistent but the new k = .1600 delta is quite a bit lighter 
(0.63 vs. 0.74). At k — .1565 we had run before only on 12 3 lattices. The new 
proton is lighter (0.84 vs. 0.89) and so is the delta (0.87 vs. 0.96). Clearly the old 
k = .1565 masses were compromised by the size of the simulation volume (as were 
the simulations at lighter valence quark mass). All am q = 0.025 baryons are also ten 
to fifteen per cent lighter than their values on 12 3 lattices. 

We do not see any oscillations in the pion effective mass at k — .1600 which we 
saw in the old doubled 16 4 running (compare Fig. ||.) There, however, the effect was 
most dramatic for the staggered valence quark pion. 

Assuming that m 2 is linear in k (as we expect from current algebra considerations), 
we can compute the critical coupling k c at which the pion becomes massless. We 
extrapolate using 

(m n a) 2 = A{- - — ) (2) 

The fit is acceptable only for the three lightest quark masses and the final numbers 
are essentially unchanged whether we use all six combinations of quarks in the pseu- 
doscalar or restrict ourselves to the three cases of degenerate quark masses. We find 
that A = 1.10(1) and k c = 0.1610(1) for am q = 0.01. These numbers are in good 
agreement with our previous results (A = 1.15(16) and k = 0.1611(1)). 

The am q = 0.025 numbers are quite different from our previous study. There we 
had A = 1.15(16) and k c = 0.1618(1). Here we have A = 1.14(1), k c = 0.1613(1) 
from the 'kind=l, WP' operator, A = 1.17(2), k c = 0.1613(1) from the 'kind=2, 
WP' operator, using only equal quark mass pions. This is such a large change that 
it cannot be due to a statistical fluctuation. When we graph the square of the pion 
mass from the old simulations (on a 12 3 lattice) and from the new simulations (on a 
16 3 lattice) we see that the new pions are consistently lighter than the old ones. The 
situation is illustrated in Fig. [|, where we also show the old and new am q = 0.01 
data. It is very strange that a finite size effect (if that is what we are seeing) would 
be stronger for heavier dynamical fermion mass. 

In Fig. |S] we present an Edinburgh plot (m p /m p vs m^/rrip). This figure also 
includes data from other simulations we have performed. Mass ratios, computed 
using correlated fits to a single exponential in each channel, are shown in Tables 
|P2| . We quantify the magnitude of hyperfine splittings in the meson and baryon sectors 
by comparing the two dimensionless quantities 

Rm = 

and 

Re = (4) 
Each of these quantities is the ratio of hyperfine splitting in a multiplet divided by 



the center of mass of the multiplet. A plot of Rm vs. Rb is shown in Fig. 10 
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Figure 8: Square of pion mass versus hopping parameter from the old am q = 0.025 
data (diamonds), new am q = 0.025 data (squares), old am q = 0.01 data (crosses) and 
new am q = 0.01 data (octagons) 
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Figure 9: Edinburgh plot for Wilson valence quarks. Data are: simulations with 
dynamical staggered fermions at (3 — 5.6 and am q = 0.01 from the 16 4 running 
(diamonds) and from the 16 3 x 32 running (crosses), (3 = 5.6, am q = 0.025 16 3 x 32 
simulations (fancy squares), and quenched simulations at (5 — 5.85 and (5 = 5.95 
(bursts). The circles show the expected infinite quark mass limit and the real- world 
point. 
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Figure 10: Comparison of baryon and meson hyperfme splitting. The two circles show 
the expected values of hyperfme splitting in the limit of infinite quark mass and from 
experiment. Points are labelled as in Fig. || 



17 



In the nonrelativistic quark model, the mass of a hadron is given by a sum of 
constituent quark masses plus a color hyperfme term, 

M h = Y,™* + ZhY<^ 1 (5) 

i ij 1 j 

where £ is twice as great for mesons as for baryons because of color [T^j From this 
model one expects the ratio Rb/Rm = 1- For all but the lightest quark mass data 
points, this is the behavior which our data shows. 

If we wish to use our spectrum results to find a lattice spacing, we can extrapolate 
particle masses to k c , fix one mass to experiment, and use this mass to infer a lattice 
spacing. We can do this for the p, proton, and A for either sea quark mass. Restricting 
our extrapolation to the lightest three valence hopping parameters, we display our 
results in Table [13|. Taking the rho as the particle whose mass is forced to its physical 
value, we have an inverse lattice spacing of 2140 MeV or 2000 MeV, proton masses 
of 1121 and 1116 MeV, and deltas at 1198 and 1302 MeV. Using the proton mass to 
set the scale, we have inverse lattice spacings of 1800 or 1685 MeV, rhos at 648 and 
650 MeV and deltas at 1008 and 1100 MeV. In all cases the proton to rho mass ratio 
is larger than experiment, and the proton-delta hyperfme splitting is too small. 



3.2 Quenched Simulations 

In an attempt to see whether any effects of dynamic fermions could be seen in the 
spectroscopy, we performed a quenched simulation with the same lattice volume as 
our dynamical simulations and with a large enough data set to overwhelm statistical 
fluctuations j| 

All fits are quite stable. We show one example of effective masses, for the (3 = 5.95, 



k = 0.1554 data set (Fig. O). The best fits to a range of points, selected using the 



histogram technique, begin at t min = 6 to 8, and are shown in Fig. 02 and Tables 14 



and ITR Mass ratios are found in Tables [16] and [T7]. Our quenched data at (3 — 5.85 



K = .1585 are consistent within statistical errors with the earlier work of Iwasaki, et. 
al-0 



While one cannot say anything about the behavior of the pion mass as a function 
of hopping parameter with two data points per (3 value, one can still extrapolate the 
square of the pion mass to zero. Doing so, we find k c = 0.1617(1), A = 1.12(4) for 
(3 = 5.85 and k c = 0.1583(1), A = 1.10(3) for (3 = 5.95. 

Finally, as a direct way of displaying any differences between our dynamical and 
quenched simulations, we show hadron masses (rho, proton, and delta) as a function of 



the pion mass in lattice units (Fig. [13]). With the possible exception of the am q = 0.01, 
k = .1600 delta, one cannot see any strong difference between spectroscopy with or 
without dynamical fermions at the parameter values used in this study. There is 
a hint in the Edinburgh plot that the nucleon to rho mass ratio in the presence of 
am q = 0.01 dynamical fermions is a bit higher than in quenched approximation. We 
display the square of the pion mass in lattice units as a function of k c — k in Fig. |14[ 
Again all the data at the lightest valence quark masses appear to lie on a (nearly) 
universal curve. 
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Figure 11: Effective masses from WP operators from quenched spectroscopy at 
(3 = 5.95, k = 0.1554. Particles in increasing order of mass are pion, rho, proton, and 
delta. 
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this paper. 
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Figure 14: Square of pion mass in lattice units for quenched and dynamical staggered 
simulations. Data are labelled as in Fig. [13] and the 16 4 
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4 Conclusions 



This concludes our program of spectroscopy for Wilson valence quarks with staggered 
sea quarks at (5 — 5.6 on 16 3 x 32 lattices. The spectroscopy we see is generally consis- 
tent with our earlier results (when performed on lattices of the same spatial volume) 
and represents an improvement over our previous results insofar as the simulation 
volume is larger. By comparing results from spatial volumes of 12 3 and 16 3 , we saw 
that for the smaller volume baryons with light valence quarks suffer from finite lattice 
size effect regardless of the dynamical fermion mass we used. This is another piece 
of evidence which suggests that sea quark properties are much less important than 
valence quark properties, in the parameter range we have studied. Perhaps we are 
also seeing the same effect on the pion mass. Of course, one can not be sure that our 
results are still not contaminated by finite volume effects on 16 3 lattices; to test that 
would require simulations in a larger volume with otherwise unchanged parameter 
values. 

When we compare our spectroscopy with dynamical fermions to quenched results 
we do not see any dramatic differences. Apparently at the parameter values of the 
simulation sea quarks simply do not affect spectroscopy above the five to ten per cent 
level. 

However, we have to say that we have always regarded this quark combination as 
something done as much for expedience as for curiosity. Wilson quarks and staggered 
quarks have very different symmetries. The next barrier in spectroscopy calculations 
will occur when lattices become large enough and quark masses become small enough 
that the rho meson mass falls below the lightest I = I, J = I -kit state on the 
lattice, or that the lightest propagating state in the rho channel is a tttc pair. These 
pions will each be made of one of the rho meson's valence quark(antiquark) and an 
ant iquark( quark) which has popped out of the sea. It would be desirable if both 
quark and antiquark have the same internal symmetry structure (for example, one 
will want to know the mass of the pion in advance), and this argues against further 
use of "hybrid" quark calculations. 
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Table 1: Fits to pseudoscalar mesons, with Wilson valence fermions and am q = 0.01 
staggered sea quarks. All fits are to "kind=l" WP propagators, unless otherwise 
indicated, and are to a single exponential. In this and following tables, numbers in 
the "kind" column for mesons refers to their quark content: 1 through 6 refer to 
hopping parameter .1320, .1410, .1525, .1565, .1585, and .1600. 
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Table 2: Fits to vector mesons, with Wilson valence fermions and am q = 0.01 stag- 
gered sea quarks. All fits are to "kind=l" WP propagators, unless otherwise indi- 
cated, and are to a single exponential. 
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Table 3: Fits to nucleons, with Wilson valence fermions and am q = 0.01 staggered 
sea quarks. All fits are to WP propagators, unless otherwise indicated, and are to a 
single exponential. 
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Table 4: Fits to deltas, with Wilson valence fermions and am q = 0.01 staggered sea 
quarks. All fits are to "kind=l" WP propagators, unless otherwise indicated, and are 
to a single exponential. 
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Table 5: Fits to pseudoscalar mesons, with Wilson valence fermions and am q = 0.025 
staggered sea quarks. All fits are to "kind=l" WP propagators, unless otherwise 
indicated, and are to a single exponential. 
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Table 6: Fits to vector mesons, with Wilson valence fermions and am q = 0.025 
staggered sea quarks. All fits are to "kind=l" WP propagators, unless otherwise 
indicated, and are to a single exponential. 
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Table 7: Fits to nucleons, with Wilson valence fermions and am q = 0.025 staggered 
sea quarks. All fits are to WP propagators, unless otherwise indicated, and are to a 
single exponential. 



28 



kind 
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Table 8: Fits to deltas, with Wilson valence fermions and am q = 0.025 staggered sea 
quarks. All fits are to "kind=l" WP propagators, unless otherwise indicated, and are 
to a single exponential. 



kind 




D m in 


D m ax 


mass ratio 


X 2 /dof 


C.L. 


(WW) 


0.1320 


12 


16 


0.997( 3) 


8.990/6 


0.174 




0.1410 
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0.985( 1) 


2.994/8 


0.935 
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11 
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6.317/8 
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0.1565 


11 
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10.740/8 


0.217 
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0.752(11) 
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0.182 
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0.558(27) 


5.007/8 


0.757 



Table 9: Fits to the ratio m^/rrip, with Wilson valence fermions and am q = 0.01 
staggered sea quarks. All fits are to "kind=l" WP propagators, unless otherwise 
indicated, and are correlated fits to a single exponential in each channel. 



kind 


K 




D m ax 


mass ratio 


X 2 /dof 


C.L. 


(WW) 


0.1320 
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16 


1.563( 8) 


11.600/16 


0.771 




0.1410 


11 


16 


1.563( 3) 


4.961/8 


0.762 




0.1525 


8 


16 


1.590( 5) 


12.880/14 


0.536 




0.1565 


7 


16 


1.591(11) 


12.520/16 


0.707 




0.1585 


8 


16 


1.597(30) 


17.370/14 


0.237 




0.1600 


6 


16 


1.573(53) 


39.600/18 


0.002 



Table 10: Fits to the ratio m^/mp, with Wilson valence fermions and am q = 0.01 
staggered sea quarks. All fits are to "kind=l" WP propagators, unless otherwise 
indicated, and are correlated fits to a single exponential in each channel. 
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kind 


^ave 


D m in 


D m ax 


mass ratio 


X 2 /doi 


C.L. 


(WW) 


0.1320 


13 


16 


1.002( 5) 


14.510/4 


0.006 


(WW) 


0.1410 


6 


16 


0.984( 2) 


25.800/18 


0.104 




0.1525 


11 


16 


0.984( 1) 


3.694/8 


0.884 




0.1565 


7 


16 


0.929( 1) 


7.176/16 


0.970 




0.1585 


7 


16 


0.845( 3) 


14.150/16 


0.588 




0.1600 


8 


16 


0.744( 7) 


14.560/14 


0.409 



Table 11: Fits to m^/rrip, with Wilson valence fermions and am q = 0.025 staggered 
sea quarks. All fits are to "kind=l" WP propagators, unless otherwise indicated, and 
are correlated fits to a single exponential in each channel. 



kind 






D m ax 


mass ratio 


X 2 /dof 


C.L. 


(WW) 


0.1320 


6 


16 


1.564( 8) 


23.060/18 


0.188 


(WW) 


0.1410 


5 


16 


1.584( 9) 


30.380/20 


0.064 




0.1525 


11 


16 


1.568( 2) 


8.207/8 


0.414 




0.1565 


11 


16 


1.596( 5) 


5.214/8 


0.734 




0.1585 


6 


16 


1.553( 9) 


21.590/18 


0.251 




0.1600 


6 


16 


1.516(15) 


21.200/18 


0.269 



Table 12: Fits to m N /m p , with Wilson valence fermions and am q = 0.025 staggered 
sea quarks. All fits are to "kind=l" WP propagators, unless otherwise indicated, and 
are correlated fits to a single exponential in each channel. 



am q 


particle 


mass 


0.01 


P 


.360(8) 




P 


.524(18) 




A 


.560(12) 


0.025 


P 


.386(5) 




P 


.558(12) 




A 


.651(17) 



Table 13: Extrapolations to k c . 
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kind 




Dmin 


D m ax 


mass 


X 2 /dof 


C.L. 


pion 


0.1585 


7 


16 


0.378( 2) 


12.382/8 


0.135 




0.1600 


7 


16 


0.273( 3) 


12.821/8 


0.118 


rho 


0.1585 


8 


16 


0.530( 6) 


2.857/7 


0.898 




0.1600 


7 


16 


0.486( 9) 


5.388/8 


0.715 


proton 


0.1585 


7 


16 


0.783(10) 


8.339/8 


0.401 




0.1600 


6 


16 


0.673( 9) 


8.113/9 


0.523 


delta 


0.1585 


8 


16 


0.852(11) 


9.302/7 


0.232 




0.1600 


8 


16 


0.757(25) 


3.628/7 


0.821 



Table 14: Fits to quenched (3 = 5.85 spectroscopy with Wilson valence fermions. All 
fits are to a single exponential. 



kind 






D m ax 


mass 


X 2 /dof 


C.L. 


pion 


0.1554 


7 


16 


0.362( 1) 


5.284/8 


0.727 




0.1567 


7 


16 


0.271( 2) 


10.564/8 


0.228 


rho 


0.1554 


5 


16 


0.486( 3) 


13.433/10 


0.200 




0.1567 


4 


16 


0.445( 5) 


19.847/11 


0.047 


proton 


0.1554 


6 


16 


0.721( 6) 


11.226/9 


0.261 




0.1567 


6 


16 


0.619( 8) 


4.814/9 


0.850 


delta 


0.1554 


6 


16 


0.777( 7) 


12.580/9 


0.183 




0.1567 


6 


16 


0.699( 9) 


7.574/9 


0.578 



Table 15: Fits to quenched (3 = 5.95 spectroscopy with Wilson valence fermions. All 
fits are to a single exponential. 



kind 


K 


D m in 


D m ax 


mass ratio 


xVdof 


C.L. 


m^/mp 


0.1585 


7 


16 


0.716( 6) 


11.050/16 


0.806 




0.1600 


7 


16 


0.573(10) 


11.740/16 


0.762 


m N /m p 


0.1585 


7 


16 


1.468(16) 


8.603/16 


0.929 




0.1600 


5 


16 


1.415(24) 


15.060/20 


0.773 



Table 16: Fits to ratios m^/rrip and m N /m p , from quenched (5 = 5.85 simulations 
with Wilson valence fermions. All fits are to "kind=l" WP propagators, and are 
correlated fits to a single exponential in each channel. 
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kind 






D ma x 


mass ratio 


xVdof 


C.L. 


m^/mp 


0.1554 


10 


16 


0.740( 8) 


5.786/10 


0.833 




0.1567 


10 


16 


0.599(16) 


7.037/10 


0.722 


m N /m p 


0.1554 


10 


16 


1.507(27) 


7.394/10 


0.688 




0.1567 


4 


16 


1.420(18) 


24.370/22 


0.328 



Table 17: Fits to ratios m n /m p and m^lm p , from quenched (5 = 5.95 simulations 
with Wilson valence fermions. All fits are to "kind=l" WP propagators, and are 
correlated fits to a single exponential in each channel. 
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